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Perturbative Matching of the NRQCD Heavy-Light Axial Current 
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A one-loop matching calculation between Lattice NRQCD and full QCD for the heavy-light axial current is 
described and the effects of renormalization on fs are discussed. 
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1. Introduction 

An important ingredient in lattice studies of B 
meson decays, is the matching between currents 
in full continuum QCD and current operators of 
the lattice theory being simulated. The GLOK 
collaboration, for instance, uses NRQCD to sim- 
ulate b-quarks together with tadpole improved 
clover light quarks [jlj. One goal is to calculate 
the pseudoscalar meson decay constant fps, 



[0\A fi \PS)=pJ PS . 



(1) 



In general, several operators of the effective the- 
ory are associated with a given full QCD current 
operator and one has, for instance, 



A, 



QCD 



I NRQCD 



(2) 



Our task is to identify the relevant operators J^' 
and calculate the matching coefficients Cj to a 
given order in a and 1/M. We present here re- 
sults of a matching calculation through 0(a/M) 
(2] and discuss their implications for fg. 



2. The Lattice Action and Current Oper- 
ators 

The lattice NRQCD heavy quark action density 
used in the perturbative calculations is 



*Talk presented by J. Shigemitsu for the GLOK Collabo- 
ration at LATTICE '97, Edinburgh, July 1997. 
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2M ' " 2M 

Cb = 1 at tree level. The one-loop contribution 
to cb is a higher order (0(a 2 /M) ) effect in the 
matching calculation and hence can be ignored 
here. For the light quarks we use the clover ac- 
tion. A few comments will be made later, on per- 
turbative calculations with Wilson light quarks. 

Heavy-light currents in full QCD have the form 
qTh. The four component Dirac spinor for the 
heavy quark, h, is related to the two component 
NRQCD heavy quark (heavy anti-quark) fields, 
ip (V), via an inverse Foldy-Wouthuysen trans- 
formation. Through 0(1/M) one has, 



U fw ^fw 



Hence at tree-level, 

j = J(°) + J(i) 
1 



— (-y-D)l 

2M V 1 n 



= qTQ 



2M 



qTh ■ D)Q 



(5) 



(0) 



Q = i(l + 7o)*FW- 

At one-loop one needs to include, for T = 7570 
and massless light quarks, a third operator 



J (2) ~(D 9 -- 7 )rQ 



(7) 
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Furthermore, on the lattice there is a discretiza- 
tion correction to 



with 

j(disc) 



r(0) 



J (0) + CU J (disc) 



(8) 



a(Dg- 7 )rg 



(9) 

CU is fixed by requiring that (J^p) projects only 
onto matrix elements that exist in the contin- 
uum theory. In perturbation theory Ca starts 
out 0(a). j( dlsc ) is the analogue of the discretiza- 
tion correction a d^P (P = pseudoscalar density) 
to the axial current in light physics emphasized 
by the Alpha collaboration. Recent studies have 
shown that this term has a non-negligible effect 
on f v around — 6.0. In heavy-light physics the 
improvement term 

j(disc) tQ j(o) also lcads t0 a 
significant correction to fs ■ 

3. Matching and One-Loop Correction 
Terms 

In order to determine the matching coefficients 
Cj we consider scattering of a heavy quark by 
the heavy- light current (qTh or the 's) into a 
light quark. The calculation is carried out in both 
full QCD and in lattice NRQCD and one requires 
that the scattering amplitudes in the two theories 
agree through 0(a/M). For the continuum QCD 
calculation we use NDR in the MS scheme. A 
gluon mass A is introduced at intermediate stages 
as an IR regulator. IR divergences cancel between 
the continuum and lattice contributions to the 
Cj's. After taking mixing among the jW's into 
account one ends up with the following one-loop 
expression, 



;a>) 



QCD 



= (l + ap )[(J^) + C A (J (disc) ) 

+ (l + apiXJ' 1 ') 

+ ap 2 (J( 2 )) (10) 

The matrix elements on the RHS are evaluated in 
the effective theory NRQCD. Writing Ca — a pd, 
one sees that there are four one-loop correction 
terms to {A ) and hence also to fps 

ap Q (J^) , «Pi(J (1) ) , ap 2 (J^) (11) 
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Figure 1. One- loop 

Corrections to a 3 / 2 f PS y/Mj^ from p (J (0) ) (Di- 
amonds) , pi (J^) (Circles) , p 2 (J^) (Crosses) 
, Pd (j( dlsc )) (Fancy Squares). Simple squares de- 
note the total one-loop correction. Vertical lines 
mark the physical B. 



and 



apd 



^j(disc)^ 



(12) 



In Fig.l we show contributions from these terms 
to a 3 / 2 fps\/Mps- Simulation results for the ma- 
trix elements were obtained on (3 = 6.0 quenched 
configurations and results are for k = n s trange 
and a — ap(l/a). One should note the relative 
importance of the discretization correction term 
a pd (j(* sc )} . At the B-meson, due to cancella- 
tion among the other one-loop corrections, this 
is the dominant one-loop term which, as we will 
see in the next section, leads to a ~ 12% reduc- 
tion of Jb relative to the tree-level result. Both 
a pd{J( dlsc ^) and a po(J^) survive into the static 
limit. The other two terms are 0(a/M) and van- 
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ish in that limit. 



The results presented sofar are for clover light 
quarks. For Wilson light quarks some modifica- 
tions are necessary. In particular, if one tries to 
calculate Ga — otpd with Wilson light quarks, 
one finds an uncancelled logarithmic IR diver- 
gence aln{a\). This divergence can only be 
removed by including clover contributions to the 
calculation. Hence, it is not possible to include 
(j( dlsc )) in a consistent way within perturbation 
theory if one is working with Wilson light quarks. 
( this is in addition to the fact that since Wilson 
light quarks have 0(a) errors one normally would 
not include 0(aa) corrections). A figure similar 
to Fig.l for Wilson quarks would only have con- 
tributions from pi(J^), i = 0,1,2. Although the 
numbers change slightly upon going from clover 
to Wilson, the qualitative feature of cancellation 
among the three one-loop corrections around the 
B meson still holds. As a result the difference 
between tree-level and one-loop corrected fs is 
smaller for Wilson than for clover. On the other 
hand there will be large scaling violations due to 
the omission of the sizeable j( dlsc ) improvement 
term. 

4. Results for fs 

Fig. 2 shows preliminary GLOK collaboration re- 
sults for, 

a 3 / 2 $ = [ a^ 2 f PS y/M^} with In(aM) in C and 
Ci set to ln(aM b ) m Zn(2.2) for all M . This en- 
sures a smooth aM — > oo limit and contact with 
HQET scaling formulas while preserving the cor- 
rect fs for the physical B meson, q* for the cou- 
pling otp{q*) is not known yet for this calculation 
and we show results for aq* = ir and aq* = 1. 
Using the p-meson mass to fix a -1 one finds, 



0.207 (9) GeV 
f Bs = { 0.189(8) GeV 
0.178 (8) GeV 



tree — level 
aq* = tt 
aq* 



= 1 



Only statistical errors are shown. After extrapo- 
lating to Kught, 

0.168 (13) GeV tree - level 
f B = { 0.153(12) GeV aq* = tt 
0.144 (11) GeV aq* = l 
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Figure 2. Heavy-Light Decay Constants at n s . 
Circles : tree-level ; Diamonds : one-loop match- 
ing with aq* — tt ; Squares : one-loop matching 
with aq* — 1. Vertical lines mark the physical B. 



In both cases, one-loop corrections are a 9% 
(14%) effect for aq* — n(l). Taking an aver- 
age, the preliminary estimate for quenched fg at 
(3 = 6.0 is, 

Jb = 0.149 (12)(t 2 5 2 )(9)GeU 

The second error comes from scale uncertainties 
and the third from higher order relativistic and 
perturbative corrections. 
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